Purity of Gaussian states: measurement schemes and time— evolution in noisy channels 
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We present a systematic study of the purity for Gaussian states of single-mode continuous variable 
systems. We prove the connection of purity to observable quantities for these states, and show 
that the joint measurement of two conjugate quadratures is necessary and sufficient to determine 
the purity at any time. The statistical reliability and the range of applicability of the proposed 
measurement scheme is tested by means of Monte Carlo simulated experiments. We then consider 
the dynamics of purity in noisy channels. We derive an evolution equation for the purity of general 
Gaussian states both in thermal and squeezed thermal baths. We show that purity is maximized 
at any given time for an initial coherent state evolving in a thermal bath, or for an initial squeezed 
state evolving in a squeezed thermal bath whose asymptotic squeezing is orthogonal to that of the 
input state. 

PACS numbers: 03.65.Yz, 42.50.-p, 03.67.Pp, 42.50.Dv 



I. INTRODUCTION 

Nonclassical features of atomic and radiation systems 
play a relevant role in quantum information, communica- 
tion and high precision measurements, as well as in many 
fundamental experiments to test quantum mechanics 
In particular, pure Gaussian states of continuous 
variable (CV) systems, such as coherent and squeezed- 
coherent states, are the key ingredients of secure optical 
communication 0, 0J|MJ6| and Heisenberg limited quan- 
tum interferometry |l H [H The characteriza- 
tion of several properties of Gaussian states has been the 
subject of intense recent work 12. 13. 14 . 15^.16,. 17-. 18]. 
stimulated by the seminal analysis on their entanglement 
properties [l3.l20|. 

Any attempt to exploit Gaussian states for quantum 
information and quantum measurement schemes must 
however face the obvious difficulty that pure states are 
unavoidably corrupted by the interaction with the envi- 
ronment. Therefore, CV Gaussian states that are avail- 
able for experiments are usually mixed states, and it 
becomes crucial to establish their degree of purity (or 
mixedness) determined by the environmental noise. In 
the present paper, we study the purity of Gaussian states 
for single-mode continuous variable systems focusing on 
two aspects: its experimental characterization, and its 
time-evolution in noisy channels. We first show that the 
joint detection of two conjugate quadratures is a neces- 
sary and sufhcient measurement to determine the purity 
of a Gaussian state with reliable experimental statistics; 
we then derive an evolution equation for the purity of 
Gaussian states in a noisy channel, considering the in- 
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stances of a thermal bath and of a squeezed thermal bath, 
and determine the evolutions that at any given time max- 
imize the purity. 

Let us refer to /i = Tr [g^] as to the purity of a 
generic quantum state g; the conjugate quantity Sl — 
{1 — fi)d/ {d — 1), where d is the dimension of the Hilbert 
space of the system under investigation, is known as lin- 
ear entropy or mixedness. In general, fi ranges from one, 
which is the value for a pure state, to fj, = 1/d for a 
completely mixed state. Throughout the paper we will 
consider CV systems, i.e. infinite dimensional Hilbert 
spaces, and therefore we will have < < 1. Since fi 
is a nonlinear function of the density matrix, it cannot 
be connected to an observable quantity if we perform re- 
peated measurements on single copies of the state. That 
is, it cannot be the expectation value of a single-system 
self-adjoint operator, nor it can be related to a single- 
system probability distribution obtained from a positive 
operator- valued measure (POVM). On the other hand, 
if collective measurements on two copies of the state 
are possible, then the purity may be measured directly 
plll23| . For instance, collective measurements of overlap 
and fidelity have been experimentally realized for qubits 
encoded into polarization states of photons |22j . 

In general, purity can be determined by the knowledge 
of the quantum state of the system, which in turn can be 
obtained by quantum tomography |24j . However, in this 
case, the statistics is usually poor, since the measurement 
of a whole quorum of observables is needed, unavoidably 
leading to large fluctuations [S^. On the other hand, if 
we focus our attention on the class of Gaussian states, 
it is indeed possible to find an operative method to ex- 
perimentally determine fi . In fact, Gaussian states are 
uniquely defined by their first two statistical moments, 
which can be measured by the joint detection of two 
conjugate quadratures, say position and momentum or 
quadrature phases of the electromagnetic field. Such a 
measurement corresponds to an estimate of the Husimi 
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Q-function Q{a) = {a\g\a), \a) being a coherent state 
of the harmonic osciUator. We will show that the mea- 
surement of the Q-function is the optimal minimal mea- 
surement for the purity, in the sense that it is necessary 
and sufhcient to determine fj. and requires the minimum 
number of observables to be measured. 

The joint measurement of two conjugate quadratures 
is possible for a single-mode radiation field as well as 
for a single atom [2fil l2^ . Remarkably, for these 
systems, the class of Gaussian states include almost all 
the states that can be reliably produced, and employed 
in communication or measurement protocols. 

Finally, we will show that the previous discussion al- 
lows to unravel the dynamics of purity only in terms of 
observable quantities. Indeed, the time-evolution of the 
purity of an initial Gaussian state in a noisy channel can 
be uniquely expressed as a function of the initial observ- 
able parameters of the input state and of the asymptotic 
observable parameters of the environment. This prop- 
erty allows then to determine and engineer optimal evo- 
lutions, i.e. evolutions that maximize the purity at any 
given time. 

The paper is structured as follows. In Section ^ we 
show how purity is related to observable quantities for 
Gaussian states, and how it can be obtained either from 
the Q-function or by three single-quadrature detection. 
In Section IIIII we present the results of a systematic 
numerical analysis, establishing the statistical reliability 
and the range of applicability of the method by means 
of Monte Carlo simulated experimental runs. We also 
show that the Q-function based determination of purity 
is a more reliable method than single-quadrature detec- 
tion. Section Hvl is devoted to derive and solve an evolu- 
tion equation for the purity of an initial Gaussian state 
in a noisy channel, both for thermal and squeezed ther- 
mal baths. We show that, even though the asymptotic 
value of purity is not related to the initial conditions, 
its behavior at finite times does depend on the initial 
squeezing and thermal excitations, and we determine the 
evolutions that maximize the purity at any finite time. 
We show in particular that purity is maximized for an 
initial coherent state evolving in a thermal bath, or for 
an initial squeezed state evolving in a squeezed thermal 
bath whose asymptotic squeezing is orthogonal to that 
of the input state. Finally, in Section we present some 
concluding remarks. 



where D{'j) = exp(7a'''— 7a) is the displacement operator, 
and Tr[0 £'(7)] is usually referred to as the characteris- 
tic function of the operator O. Let Oi and O2 be oper- 
ators that admit regular Wigner representations Oi(a) 
and 02(a)- Then the trace Tr [Oi O2] can be computed 
as an integral over phase space according to 



Tr[Oi O2] = TT / d^a Oi(a) 02(a) 
Jc 



(2) 



From now on, we will move to the phase-space vari- 
ables X and p, corresponding to quadrature phases x = 
(a + a^)/^/2 and p = i(a^ — a)/V2 of the field a, whose 
expectation values (x) = x and (p) = p are related to a 
by a = (a; -|- ip) / \/2. 

The Wigner representation W{a) of the density matrix 
p of a quantum state is referred to as the Wigner function 
of the state. The class of Gaussian states is defined as the 
class of states with Gaussian Wigner function, namely 



Wix,p) 



(3) 



where X is the displaced vector X = (.t — xo,p — po) and 
cr is the covariance matrix 



where xi — x,X2 — P- The density matrix of the most 
general Gaussian state can be written as j30i | 



g ^ D{ao)S{r, Lp)vnS^{r, Lp)D^{ao) , 



(4) 



where — (xo + «po)/V2, and Vfi is a thermal state with 
average photon number n: 



|fc)(fc|. 



1 + n 



St 



fc=0 



Z3(ao) denotes the displacement operator and S{r,if) — 
exp(ir e~*^'^a^ — e^'^'^a^'^) the squeezing operator. A 
convenient parametrization of Gaussian states can be 
achieved replacing the CTi^'s by n, r, ip, which have a 
more direct phenomenological interpretation. By apply- 
ing the phase-space representation of squeezing |29l l3^ , 
the following relations are easily derived 



II. PURITY OF GAUSSIAN STATES 

We begin by reviewing some fundamental properties of 
the Wigner phase-space representation which will be 
useful throughout the paper. The Wigner representation 
of an arbitrary operator O is defined as follows 
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cosh(2r) — sinh(2r) cos(2iy9) 
cosh(2r) + sinh(2r) cos(2iy9) 
sinh(2r) sm{2(p) . 



Exploiting Eq. one can write 

M = Tr[e2]^ J f [ dxApW\x,p), 



(5) 



(6) 
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so that, for a Gaussian state 
1 



2^Det[c 



(7) 



In terms of n, r and ip, Eq. ||7J) can then be recast as 

ISlIsi 



1 



2n+ 1 



(8) 



Eq. © shows that the purity of a generic Gaussian state 
depends only on the average number of thermal photons, 
as one should expect since displacement and squeezing 
are unitary operations. Therefore, the measurement of 
the purity of a Gaussian state is equivalent to the mea- 
surement of its average number of thermal photons. 

As the last step in connecting fi to observables we re- 
port the expression of the aij 's in terms of the Q-function 
Q{a). This follows from the antinormally ordered expres- 
sion of the second moments. We have, for instance 



which, in terms of phase-space variables, corresponds to 



„2 1 



2 . Therefore, we eventually get 
(x ) = Ti-[q ] 



/ [ dxdpQ{x,p) (x^ - i) , 



where we have moved from variables a and a to variables 
X and p, previously defined. In much the same way, we 
obtain 



dxdpQ{x,p){p^ -^), (9) 



-(xp + px) — dx dp Q{x,p) X p . 

2 Jr Jr 



(10) 



Since first moments are naturally antinormally ordered, 
evaluation of first moments of quadratures is easily ob- 
tained, and the cry 's can be eventually computed. 

Gaussian states may be effectively characterized as well 
by single-quadrature measurements obtained by balanced 
homodyne detection [s^l ■ Thus a question arises whether 
or not one really needs to resort to joint measurement of 
two conjugate quadratures to determine the purity. In 
particular, since Gaussian states are fully characterized 
by the first and second moments, it suffices to measure 
the rotated quadrature xg = (a''' e'^-|-a e~*^)/-\/2 for three 
different values of 9 to have a complete characterization 
of the state, including the measure of its purity. This fact 
can be proven by reminding that the probability distri- 
bution p{x, 9) of a measurement of xg on a state of the 
form Q is a Gaussian centered in xo = Re[ao e~*^], with 
variance 
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e-2-cos2( 



^) + e^"- sin\9 - if)] . (11) 
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FIG. 1: Effect of the number of data on the Q-function based 
determination of purity for Gaussian states: results from 
Monte Carlo simulated experiments. On the left graph we 
plot the determination of the purity ^ versus the number of 
data Nx for a squeezed thermal state with parameters given 
by ao = 0, ^ = 0, r = 1.5, corresponding to sinh^ 1.5 ~ 4.5 
mean squeezed photons, and a mean number of thermal pho- 
tons n — 0.5. Black circles are the estimated values of purity 
based on the Q-distributed statistics, vertical bars are the 
experimental errors (confidence intervals); for a large number 
of experimental data the errors quickly fall well within the 
black circles of the estimated values. The theoretical value of 
purity for all the simulated experimental runs is fi = 0.5. On 
the right graph we report the relative errors fJ. versus the 
number of data for the same squeezed thermal state. 



By measuring three quadratures we directly obtain the 
purity /i by comparison of variances. By choosing 9 = 
0, 7r/2, 7r/4 we have 

M = [4cr7r/4(CTo + - - (cto " ' • (12) 

In the next Section, we will compare the two different 
experimental schemes on the basis of Monte Carlo simu- 
lated experiments. 



III. MONTE CARLO SIMULATED 
EXPERIMENTS 

As we have seen, in order to evaluate the tr^j's and 
then the purity, we need to estimate averages over the 
Q-function. These estimates can be obtained if one dis- 
poses of data distributed according to the Q-function 
Q(x,p) itself. Indeed, such a distribution can be experi- 
mentally reconstructed for a single-mode radiation field 
through heterodyne js^ , eight-port homodyne [s^ [s^ 
or six-port homodyne detectors jsl], and for atoms by 
coupling the atom with two light fields and measuring 
the corresponding phase-shifts . 

In order to test the effectiveness of the proposed 
scheme, we have performed a systematic numerical analy- 
sis by means of Monte Carlo simulated experiments. The 
simulations are needed to show the actual independence 
of the method on the squeezing and displacing parame- 
ters, in compliance with Eq. ijHJi. Moreover, they provide 
a crucial test on the actual possibility of getting reliable 
(i.e. with reduced fluctuations) determinations of /i in 
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FIG. 2: Effect of the number of data on the determina- 
tion of purity for Gaussian states by single-quadrature de- 
tection: results from Monte Carlo simulated experiments of 
three quadratures xq, x„/2, x^^4. On the left graph we plot 
the determination of the purity fi versus the number of data 
Nx for the same squeezed thermal state of Fig.0 Notice that, 
in this instance, the total number of data corresponds to 
Nx/S detections for each quadrature. Black circles are the 
estimated values of purity based on the balanced homodyne 
statistics, vertical bars are the experimental errors (confidence 
intervals); for a large number of experimental data the errors 
quickly fall well within the black circles of the estimated val- 
ues. The theoretical value of purity for all the simulated ex- 
perimental runs is fi = 0.5. On the right graph we report the 
relative errors Ajj,/ /j. versus the number of data for the same 
squeezed thermal state. 



realistic experimental settings and even for most unfa- 
vorable states. 

The purity /i and its dispersion A/i have been evalu- 
ated from samples of the Q-function, varying the values 
of the parameters of the simulated Gaussian state. Be- 
sides n, r, if, and ao, the experimental determination of 
the Q-function depends on the number of collected 
data. 

We find that /i and A/i are essentially independent 
on the complex displacement parameter ao and on the 
squeezing angle ip. On the other hand, Afi does depend 
on n and r, decreasing with increasing n and increasing 
with increasing r. 

In Fig. n we report the determination of purity for 
a strongly squeezed thermal state as a function of the 
number of data. The error on purity is of the order of a 
few percent for samples made of — 10^ data. 

In order to compare the determination of /x by the 
Q-function with that coming from single-quadrature de- 
tection, we have simulated the measurement of three 
quadratures xg, 9 — 0, 7r/2,7r/4 by balanced homodyne 
detection. In Fig. [3 we report the estimated purity [using 
Eq. (|12|l ] for the same strongly squeezed thermal state of 
Fig. n as a function of the number of data. Some fea- 
tures are immediately evident. First of all one can see 
that the determination is biased: in the present case the 
estimated ji is always larger than the true value, while 
the opposite case occurs by inverting the phase of the 
squeezing. Therefore the method is very sensitive to the 
choice of the phase. Moreover, the relative error is not a 



FIG. 3: Effect of squeezing on the Q-function based deter- 
mination of purity for Gaussian states: results from Monte 
Carlo simulated experiments. On the left graph we plot the 
determination of purity versus the squeezing parameter r for 
Gaussian states with the other parameters fixed at a = 0, 
ip = Q, and n = 0.5. Black circles are the determined values 
of purity based on the Q-distributed statistics, vertical bars 
are the experimental errors (confidence intervals). For small 
r the errors are within the black circles. The theoretical value 
of purity for all the states is ^ = 0.5. On the right graph we 
report the relative errors Afi/fj, versus the squeezing param- 
eter for the same set of experiments. The number of data in 
all simulated experiments is Nx — 3 • 10^. 



smooth function of the number of data i.e. the method 
is not statistically reliable as the joint-measurement one. 
This is again due to the remarkable dependence of the 
variances on the phase of the squeezing, a dependence 
which is instead smoothed out in the measurement of 
the Q~function. Summing up, for some specific states 
(as the example considered here) single-quadrature de- 
tection may be asymptotically even more efficient than 
the heterodyne one. However, in general, the number of 
data needed for the relative error to be below the joint- 
measurement level is strongly state-dependent. We con- 
clude that the measurement of the Q~function is statisti- 
cally more reliable and thus more suited for a systematic 
analysis of the purity of Gaussian states. 

Let us now go back to the analysis of the Q-function 
determination of purity. A smaller number of data is 
needed to obtain a given precision for states with smaller 
squeezing. The effect of the squeezing parameter on the 
determination of purity is illustrated in Fig.|21 where we 
report n and Afj,/^ versus r for Gaussian states with a = 
0, (/? = 0, and n = 0.5, and for a number of data — 
3 • 10''. Notice that in Fig. O the range of r corresponds 
to a quite large number of mean squeezed photons < 
sinh^ r < 15. 

In the deep quantum regime, i.e. for small n, fiuctu- 
ations of fi become more relevant. This is not surpris- 
ing, since /i is a highly nonlinear function of the second- 
order moments. However, simulations show that even for 
highly squeezed (up to ~ 15 mean squeezed photons) and 
slightly mixed (down to n ~ 0.1) states, realistic experi- 
mental conditions allow a statistically reliable determina- 
tion of fi that complies with the theoretical expectation 
©, up to an error of a few percent. In Fig. ^ we plot 
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FIG. 4: Effect of thermal photons on the Q-function based de- 
termination of purity for Gaussian states: results from Monte 
Carlo simulated experiments. On the left graph we plot the 
determination of purity versus the value of n for Gaussian 
states with the other parameters fixed at a = 0, = 0, and 
r = 1.0. Black circles are the determined values of purity ac- 
cording to the Q-function statistics, and vertical bars denote 
the experimental errors (confidence intervals); the latter are 
within the black circles for essentially all values of n. The 
solid line reports the theoretical values of /x. On the right 
graph we report the relative errors Afi/fj, versus n for the 
same set of experiments. The number of data in all simulated 
experiments is Nx — 10''. 



the determination of purity for different squeezed ther- 
mal states as a function of the average number of thermal 
photons n, for samples made of — 10^ data. 

From the above analysis we conclude that the joint 
measurement of two conjugate quadratures provides a 
statistically reliable method to determine the purity of a 
generic Gaussian state. This is best achieved with exper- 
imental schemes that involve data distributed according 
to the Husimi Q-function, such as heterodyne and multi- 
port homodyne detection schemes. 



IV. EVOLUTION OF PURITY IN A NOISY 
CHANNEL 

Let us consider the time evolution of an initial, pure 
or mixed, generic single-mode Gaussian state in presence 
of noise and damping (and/or pumping) toward a final 
squeezed thermal state. If F^^ is the photon lifetime in 
the noisy channel, the evolution of a state is described, in 
the interaction picture, by the following master equation 

g = ^NL[a^g+^{N+l)L[a]g 



M D[a]g + M D[a'']g 



(13) 



where the dot stands for time-derivative and the Lind- 
blad superoperators are defined by 

L[0]g = 20gO^ -O^Og- gO^'O, (14) 



M is the correlation function of the bath (which is usually 
referred to as the squeezing of the bath) ; it is in general a 
complex number M = AIi+iM2, and M denotes its com- 
plex conjugate, while is a phenomenological parameter 
related, as we shall see, to the purity of the asymptotic 
state. Positivity of the density matrix imposes the im- 
portant constraint |Mp < N{N + 1). At thermal equilib- 
rium, i.e. for M = 0, N coincides with the average num- 
ber of thermal photons in the bath. The master equation 
(|13|1 can be transformed into a Fokker-Planck equation 
for the Wigner function W{x,p,t). U sing th e differential 
representation of the superoperators |33,|33 in Eq. I|13|l . 
the corresponding Fokker-Planck equation reads as fol- 
lows 



xdx + pdp 



2N + 1 



(a; 



,2 

XX 



ppj 



+ Ml (C - dip) + 2M2d,pj W{x,p, t).(16) 

For a general single-mode Gaussian state of the form Q 
one has, in compact notation 




+ MiSxA^) + ^l2SxA«)jW, (17) 

where Sxa-i'j) denotes the seralian (or siralian) operator, 
a scalar function of the matrix 7 given by Sxa-i'j) = 
X(T-^jcr-^X'^ - Tr[7cr-i]. The displaced vector X and 
the covariance matrix cr have been previously defined, 
whereas I, A, B form a basis in the space of 2 x 2 real 
symmetric matrices: 



1 

1 



For any given real matrix 7 and generic Gaussian states, 
the seralian operator shows the remarkable property: 















A = 





dx / dpSx^h)W{x,p,t)^Q. 



(18) 



It can be easily shown that this property assures that the 
last three terms of Eq. H17|) [corresponding to diffusion 
terms in the Fokker-Planck equation 1)16(1 ] do not enter 
in the time-evolution equations for the first statistical 
moments xq and pq. Such evolution is governed by the 
drift terms and is described by the following equation for 
the vector Xq = 



dx / dp 



x\ W 
P 



-Xn 



(19) 



D[0]g = 20gO -OOg- gOO 



(15) 



First moments are damped through the noisy channel: 
this effect should be expected since it is the mathemat- 
ical evidence of the absorption of the state's coherent 
photons. 
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The evolution of the covariance matrix of the state can 
be described by monitoring different sets of variables. A 
good choice of variables is given by the a^-'s, in terms 
of which the evolution equations decouple. The relations 
between the variables aij and the variables fi, r, and ip are 
given in Eqns. ^ and (0. Here we recall some further 
relations that will be useful in the following: 



Det[<T] = axx(^pp — CT, 



(2n+ 1) 



xp 



Tr[cr] = axx + CTpp = (2?! + 1) cosh(2r) = 



1 

cosh(2r) 



sinh(2r) cos(2iy9) 



(20) 



(21) 



(22) 



As we have seen, in the Wigner phase-space picture the 
expectation values can be computed as phase-space in- 
tegrals. The first-order evolution equation for the co- 
variance matrix <r is thus obtained by straightforward 
integration, and reads: 



with 



r {(Too - a-) 



(2jV+l)+2Mi 



(2jV+l)-2Mi 



(23) 



(24) 



In fact, integration of Eq. (|23|) yields 



ait) 



(1 



(T(O)e 



-rt 



(25) 



Eq. H25(l shows a simple example of a Gaussian com- 
pletely positive map [i^. The Gaussian character of the 
evoluted Wigner function can be proven, a posteriori, by 
verifying that a function of the form Q , with covariance 
matrix given by Eq. H25(l , indeed solves Eq. . In order 
to be a bona fide covariance matrix, cr{t) must satisfy the 
usual condition encoding the x — p uncertainty relations 

MM 



T(i) + -J>0, 



with J = 



1 

-1 



(26) 



It is promptly seen that such a condition is satisfied 
at any time by the convex combination giving (j{t) in 
Eq. (|25|1 iff (Too is a legitimate covariance matrix. This 
last requirement is assured by the necessary constraint 
N{N -f 1) > |Mp that guarantees positivity of the den- 
sity matrix. 

By introducing 



Mo 



(2iV+ 1)2 -4|M|^ 



'1/2 



and exploiting Eqns. (|20H22|) we can eventually express 
/i, r and as functions of time 



The matrix crao, determined by the bath parameters 
alone, turns out to be the asymptotic covariance matrix. 



IJ.{t) = /Xo 



4(1 



-2Tt 



2mo( 



Vl + 4M?o|M|2cosh(2ro) 



H-2sinh(2ro)(MiCOs(2(po) - M2 sin(2v5o)) ) (l 



-rn „-rt 



-1/2 



cosh[2.(0] = N'^^^UM?i'~^-'l^,-r.^J^M^ 



ia.n[2ip{t)] = 



Moo Mo 

M22^o (1 - e-i'*) + sinh(2ro) sin(2.^o)e"" 
-Afi2/io (1 - e-r*) + sinh(2ro) cos(2(^o)e" 



-rt 



(27) 
(28) 
(29) 



where mo, '^o and ip^ are, respectively, the initial purity and the initial squeezing parameters. 

Let us first consider the case M = 0, for which the initial state is damped toward a thermal state with mean photon 
number j2£l, 42] . In this case, see Eq. H29|l . ip is constant in time and does not enter in the expression of m- The 
corresponding solutions for /i(t) and r{t) read then as follows 

2 n-i/2 

\^ — I -r ^ 

Moo 



Mo 



M(t) 

cosh[2r(t)] = ij.{t) 



-^t^- , o A^o „_rt „-rn cosh(2ro) + e'^^* 



2^(l-c-^^Y + 2^^-^' (1-e 
MSo / „ V 



1 - c 



-r* 



Moo 



_P( cosh(2ro) 
Mo 



(30) 
(31) 
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FIG. 5: Plot of the purity /i for an initially pure Gaussian 
state {fio = 1) in different non-squeezed (M — 0) noisy chan- 
nels, evaluated at time t = F"^, as a function of the initial 
squeezing parameter ro. From top to bottom, the value of 
the mean thermal photon number TV that characterizes the 
different channels isA'^^0, A'^ = 0.5, and A*' = 1 respectively. 



The quantities n{t) and r{t) in Eqns. H30|) and (|31|l solve 
the following system of coupled equations 



r 



cosh(2r) 



r = — sinh(2r) , 

2 /ioO 



(32) 



which, in turn, can be directly found working out the 
basic evolution equation /i = 2TT[g g] as a phase-space 
integral and exploiting Eqns. H2(JH22|I . It is easy to see 
that, as i oo, fi{t) ^ fj,^ = (27V + 1)-^ and r{t) 0, 
as one expects, since the channel damps (pumps) the 
initial state to a thermal state with mean photon number 
iV. Therefore, the only constant solution of Eq. (|32|) is 



M = Mo 



0, i.e. only initial non-squeezed states are 



left unchanged by the evolution in the noisy channel. In 
fact Eq. H3U|) shows that /i(t) is a decreasing function 
of vq: in a non-squeezed channel (M = 0), a squeezed 
state decoheres more rapidly than a non-squeezed one 
(see Figs. El and El). Let us consider, for instance, an 
initially pure state in a channel with N = 1 (so that 
Moo = |); after a time t = F^^, the ratio of the purity of 
a state with ro — 1.5 to the purity of a state with tq = 
is 53.7%. This dependence could therefore be relevant for 
practical purposes. The optimal evolution for the purity, 
obtained letting r = in Eq. (|30|l . reads 



Mo Mo 



Mo 



(Mo 



Moj 



(33) 



Obviously, ii{t) is not necessarily a decreasing function 
of time: if mo < Moo then the initial state will undergo 
a certain amount of purification, asymptotically reaching 
the value Moo which characterizes the channel, as shown 
in Fig. El In addition, fi{t) is not a monotonic function 
for any choice of the initial conditions. Letting jj, = in 



Eq. H32() . and exploiting Eqns. (|30|l and H31|) . one finds 
the following condition for the appearance of a zero of 
/i at finite positive times: cosh(2ro) > max(^,^). 
If this condition is satisfied, then /i(i) shows a local ex- 
tremum, in fact a minimum since, differentiating the first 
of Eqns. and letting /i = 0, one obtains jj. > 0. This 
behavior is shown in Fig. El 

Let us now treat the more general instance M ^ Q 
of a squeezed thermal bath. Recalling the definition of 
/ioo and exploiting Eqns. H27l - I29|l . one easily finds the 
asymptotic values of the physical parameters Mi i^ and ip 



1 



v/(2iV+l)2-4|Af|2 



sh(2roo) = V1 + 4mLI*^P, 



tan(2(^oo) = 



(34) 



(35) 



(36) 



These values characterize the squeezed channel. Eq. (|34f) 
shows that, if M 0, then N is not simply the mean 
thermal photon number n of the asymptotic state. One 
has: 



1)^ 



-4|M|2- 1 



N 



In order to understand the dynamics of purity when 
M ^ 0, it is convenient to write again the expression 
(|27|l for m(0; using Eqns. H35() and 136() to switch from 
the complex parameter AI = Mi + iM2 to the asymp- 
totic values of the squeezing parameters Too and (poo one 
obtains 




FIG. 6: The purity jj, for various Gaussian states evolving in a 
channel with N — 0.5, M = 0, as a function of time. Time is 
dimensionless and measured in units of F"^ . The upper curve 
refers to an initial pure coherent state (ro = 0, fj,o ~ 1), the 
central curve to an initial pure squeezed vacuum (ro = 1.5, 
jio ~ i), and the lower curve to an initial thermal state with 
ro = and jj.o — 0.05, i.e. no = 9.5. 
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Mo 



Mo 



,-2rt 



+ 2-^ ( cosh(2roo) cosh(2ro) + sinh(2roo) sinh(2ro)( cos(2(^oo — 2(^o))) 



-1/2 



(37) 



We see from Eq. (|l-i7ll that is a monotonically decreasing function of the factor cos(2(y9oo — 2(po), which gives the 
only dependence on the initial phase (fo of the squeezing. Thus, for any given (poo characterizing the squeezing of the 
bath, ifQ = ipao + § is the most favorable value of the initial angle of squeezing, i.e. the one which allows the maximum 
purity at a given time. For such a choice, ^It) reduces to 



Mo 



1-e 



-rt\ 



+ c 



-2rt 



2-^ cosh(2ro, - 2ro) (l - e"") e"^* 

Mo 



-1-1/2 



(38) 



This is a decreasing function of the factor cosh(2roo — 
2ro), so that the maximum value of the purity at a given 
time is achieved for the choice rg — , and the evolution 
of the purity of a squeezed state in a squeezed channel is 
identical to the evolution of the purity of a non-squeezed 
state in a non-squeezed channel expressed by Eq. H33|) 
and illustrated in Fig. |H1 

In conclusion, for the most general instance of a chan- 
nel characterized by arbitrary fioo, r^o, foo and F, the 
initial Gaussian state for which purity is best preserved 
in time must have a squeezing parameter = Too and a 
squeezing angle ipQ = <y5oo + f , i.e. it must be antisqueezed 
(orthogonally squeezed) with respect to the bath. The 
net effect for the evolution of the purity is that the two 
orthogonal squeezings of the initial state and of the bath 
cancel each other exactly, thus reproducing the optimal 
purity evolution of an initial non-squeezed coherent state 
in a non-squeezed thermal bath. 

V. CONCLUSIONS 

We have shown that the purity of Gaussian states for 
continuous variable systems can be operationally deter- 
mined by the joint measurement of two conjugate quadra- 
tures. In order to perform such a measurement, the mini- 
mal, necessary and sufficient requirement is that the mea- 
surement apparatus records data distributed according 
to the Husimi quasi-probability function. We have then 
verified by Monte Carlo simulated experiments the statis- 



tical reliability of the associated measurement schemes, 
thus proving the possibility of an experimentally real- 
izable characterization of the purity of Gaussian states. 
We have compared as well the scheme based on the Q- 
function with the one based on single-quadrature detec- 
tion, and showed that the former provides a more reliable 
statistics. Moreover, we have derived an evolution equa- 
tion for the purity of Gaussian states in noisy channels, 
both in the case of a thermal and of a squeezed thermal 
bath. Our analysis shows that the purity is maximized at 
any given time for an initial coherent state evolving in a 
thermal bath, or for an initial squeezed state evolving in 
a squeezed thermal bath whose squeezing is orthogonal 
to that of the input state. We have focused our attention 
on the purity of single-mode Gaussian states. The time- 
evolution of the purity for specific initial non Gaussian 
states of great physical relevance can be studied, as well 
as the extension to Gaussian states of multimode sys- 
tems, both pure and mixed. These topics are currently 
being explored and will be the subject of forthcoming 
work. 
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